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Abstract 

The Observed Perlmutter-Riess Acceleration (OPRA) implies that 
the expansion of the Universe is currently increasing and is motivation 
for the so-called phantom energy models. We consider the dynamics of 
phantom scalar field models. An important physical time constraint, 
which can be used to rule out many cosmological models, is obtained 
from the condition that all forms of energy density, including the field 
causing OPRA (e.g., the phantom field), must be non-negligible for an 
extended period, which is conservatively estimated to be of the order 
of a few Gyr. We find that this physical time constraint cannot be 
satisfied in conventional phantom cosmological models. 


1 ACT I: OPRA 

Evidence from supernovae observations .TJ strongly suggest the possibility of 
late-time accelerated expansion of the Universe; the Observed Perlmutter- 
Riess Acceleration (OPRA). The energy density is usually described by an 
effective “equation-of-state” parameter r y e ff~ 1 = p/p, the ratio of the energy 
pressure p to its energy density p. A value 7 e // < 2/3 is required for cos¬ 
mic acceleration. The most recent observations, including supernovae, cos¬ 
mic microwave background experiments, quasar-lensing statistics and galaxy 
clustering observations (see 0 and references therein), taken together sug¬ 
gest a Universe with 7 e ff satisfying 7 e ff < 0; this is consistent with 7 e ff = 0, 
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but it has been argued that with certain prior assumptions the best fit is 
actually for 7 e ff < 0 0 EI. It is therefore important to look for theoretical 
possibilities to describe dark energy with 7 e ff < 0 . 

A wide class of scalar held cosmologies have been utilized to model dark 
energy. Quintessence models [ 0 ] lead to 0 < 7 e // <2/3. The simplest expla¬ 
nation for dark energy is a cosmological constant, for which 7 e ff = 0. Matter 
with 7 e ff < 0 has been dubbed “phantom energy” j 2 ], and has received in¬ 
creased attention recently 0. Specific models with non-minimally coupled 
scalar fields may lead to phantom energy; perhaps the simplest alternative 
is provided by a phantom scalar held with negative kinetic energy. Theorists 
have also discussed stringy phantom energy |BJ and brane-world phantom 
energy jjj. 

Models with a constant 7 e ff < 0 lead to a future singularity commonly 
called the Big Rip. This singularity is typically characterized by a divergent 
pressure and acceleration in a finite proper time (although there are rip-free 
phantom cosmologies and there are also sudden future singularities that are 
less severe j§]). The energy density of phantom held increases and eventually 
violates the dominant energy condition. The physical consequences in this 
scenario have been discussed in [3]. Expanding universes that come to a 
violent end after a finite proper time at a future Big Rip singularity have a 
number of additional undesirable properties |5j. 

In this paper we take a conservative view and we will only analyze the 
cosmological dynamics of the phantom held and relate it to observations. 
The dynamics of phantom cosmology has been discussed previously urn- 

2 ACT II: The Phantom 

The Einstein equations in a Robertson-Walker spacetime containing phantom 
matter and a separately conserved perfect fluid are 


2k 

— (Pb + Pb + Pph + Pph ) + 7^2 

( 1 ) 

9 / > 3 k 

3 H = (pb + Pph )- 7 , 

GT 

( 2 ) 

pb + 3 H(p b +p b ) = 0, 

( 3 ) 


where a(t) is the expansion scale factor, k is the curvature parameter, and 
H = a/a is the Hubble expansion rate. In the above we have assumed that 
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the perfect fluid satisfies the barotropic equation of state 

Pb = (7 - 1 )pb, (4) 

where 7 is a constant which satisfies 0 < 7 < 2 . 

In the best motivated models the phantom cosmologies are generated by 
a (separately conserved) scalar held with negative kinetic term. The energy 
density and pressure of the held are given by 


Pph — 4 * ^( 0 ), 

(5) 

p P h = 

(6) 


The corresponding equation of state parameter is now given by ( 7 p h — 1) = 
P<!>/Pb < 0, f° r P<f> > 0 - The effective Klcin-Gordon equation reads 

dV 

o - 3//o - — 0. (7) 

We shall consider a self-interacting scalar held with an exponential potential 

4" = V 0 e**, ( 8 ) 

where Vo and k are positive constants. We note that a qualitative analysis 
of the standard spatially homogeneous and isotropic scalar held models has 
shown that there exist scaling solutions which can act as unique late-time 
attractors in these models m- 


2.1 Dynamics 

Defining 

^ 0 T V 

y/GH ’ 3 H 2 

and the new logarithmic time variable r by 


Q b 


Pb 

3 H 2 ' 


(9) 


( 10 ) 
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the governing equations can be written as the three-dimensional autonomous 
system: 




T' 


9 ' 



9b — jj (2 + 2 <F 2 + T) , 
9b — ^(24> 2 + 'h — 1) , 
-^(24> 2 + T) , 


( 11 ) 

( 12 ) 

(13) 


where a prime denotes differentiation with respect to r, and the Friedmann 
equation becomes 

1 - 9 b - (T - 4> 2 ) = - k/a 2 H 2 (14) 


The physical region of the state-space is constrained by the requirement 
that the energy densities be non-negative; i.e., > 0 and 9 p h = (T — <f> 2 ) > 

0. Furthermore, from equation Cl we hnd that the variables 9b and 9 p h 
are bounded for k = 0 and k = — 1, but not for k = +1. The individual 
variables need not be bounded. Note that 9 p h = 0 is not an invariant set. 

Geometrically the zero-curvature models (k = 0) are represented by 9b — 
<f> 2 + T = 1, in the (<f>, T, f2 b ) state-space. Defining /C = 1 — 9 b — (T — <f> 2 ), 
we see that 


K! 


1C 



2(2<f> 2 + T) , 


(15) 


so that /C = 0 is an invariant set. 


2.1.1 Equilibrium points at finite values 

The equilibrium points ($ 0 , T 0 , 9 b o) in the physical phase space at finite 
values are: 

M (0,0,0): eigenvalues [—2, 2, (2 — 3y)] (saddle). The exact solution corre¬ 
sponding to this equilibrium point is the Milne model. 

F (0, 0,1): eigenvalues [—1(2 — 7 )), 3y, (3q — 2)] (saddle). The exact solu¬ 
tion corresponding to this equilibrium point is the flat Friedmann model. 

R 


^,1 + yj-jO): eigenvalues 
The exact solution corresponds 


-(2 + K 2 ), -3(1 + ^), -(3y + ft 2 ) (sink). 
;o a future late-time singularity called the 


Big Rip singularity. In this power-law solution 

a = a 0 (—t)~ 2 ^ 2 , 


(16) 
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V 


2(6 + k 2 ) ^ 

K 4 


(17) 


-2 

0 = — In \t\. (18) 

This solution represents a universe in which the scale factor grows mono- 
tonically till a future Big Rip is reached at t — 0. In our models, p > 0 
but the weak energy condition is violated because p + p ~ —0 2 < 0. From 
flEJ), the effective barotropic index for the power-law phantom held solution 
is 7 p h = —4> 2 /Pph = — ft 2 /3 which becomes negative while both the energy 
density p p h and the pressure p p h of the held diverge (at t = 0 ). 

We note that there are no equilibrium points corresponding to matter 
scaling solutions, in which both ordinary/dark matter and phantom matter 
are non-negligible ($ 0 2 + T 0 ) ^ 0,h2b 0 7 ^ 0 [Hi]. This was noted in EH- In 
addition, it was shown El that there are no matter scaling solutions in two 
scalar held models (with exponential potentials). 


2.1.2 Equilibrium points at infinity 

We need to consider what happens if the variables diverge. This can occur if 
H = 0, thus leading to bouncing models [13]. Let us therefore assume that 
k < 0 and $, T —> oo. From the equations of motion we note that this can 
only happen to the past. In fact, the variables (including Vt b ) will in general 
diverge to the past; asymptotically, we find 

<F 2 , T, fl b oc —-—, r > r 0 , (19) 

T~ To 

and hence, the variables blow up in finite r. This blow-up is only an artifact 
of the variables chosen and results because the universe experiences a bounce. 
In principle, alternative dimensionless variables can be chosen to render the 
state space compact, even when the universe bounces. More details are 
presented in [ITT)] . 

A more serious problem with these models is the behaviour of the phan¬ 
tom energy density. We find that Vt p h = T — $ 2 goes negative into the past, 
and that the models, without modifications, are not physical (HI- 


2.2 Alternative Phantoms 

There is an alternative approach to obtaining Big Rip singularities in which 
the phantom energy is modelled as a (separately conserved) perfect fluid with 
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effective equation of state p p h = (x P h ~ 1 )p P h, where 7 p h < 0 is a constant. 
Defining normalized variables, all equilibrium points are found to have either 
Dfe or fl p h zero, and hence again there are 110 matter scaling solutions. Indeed, 
defining x = (^ph~^b)/(^ P h+^b), we find that x is monotonically increasing 
from y = —1 to x — +1. These models cannot be physical. 

It has been noted that matter scaling solutions are not possible in phan¬ 
tom scalar field models with other potential HU. However, scaling solutions 
are possible if there are interactions allowed between the phantom field and 
the background matter fields HU; the physical consequences of these inter¬ 
acting models remain to be fully investigated. 

3 ACT III: the Phantom Unmasked? 

The physical models need to satisfy a number of observational constraints. 
In general, the Universe is approximately flat (but has non-zero curvature), 
the Hubble time Ho 1 corresponds to h = 0.71 ± 0.06 d. and it is strongly 
indicated that presently 5 percent of matter is due to baryons/ordinary mat¬ 
ter, 25 percent dark matter, and 70 percent dark energy. The coincidence 
problem (or fine tuning problem) asks why are the matter density, scalar 
density, and phantom density comparable at the present epoch. 

In particular, the matter and scalar fields must be non-negligible for an 
extended period in order to be consistent with galaxy formation, the existence 
of dark matter and OPRA. This is difficult to achieve in regular cosmology, 
although it can be accomplished in matter plus scalar field models due to the 
existence of a matter scaling solution PI, which corresponds to an equilib¬ 
rium point in a dynamical systems analysis (but not necessarily an attractor; 
usually in the presence of curvature this equilibrium point is a saddle, and if 
trajectories move close to this saddle, the trajectories will stay close for an 
arbitrarily long period of time before subsequently evolving away). 

The transition from deceleration to acceleration occurs at redshift zt > 
0.4 HUGH, and can be higher depending on the parameters of the model; this 
redshift corresponds to t T = 0.6t 0 , where t 0 is the present time. Consequently, 
the time period during which the phantom field must be uon-negligible, t p h, 
satisfies t p h > to — tr ~ \{H o _1 ). 

I 11 the theory of structure formation the present day galaxies and galaxy 
clusters formed due to gravitational instability of initially small primordial 
density fluctuations PI- Thus the time U for which the matter is non- 
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Figure 1: The evolution of the phantom energy density, f } p h , and the energy 
density of matter, in terms of the logarithmic time r. 


negligible must be at least as long as the age of globular clusters ESI : con¬ 
sequently, t b > 12.0 Gyr. We also note that t ( i m > to ~ tree- 

We shall not impose all of these constraints here, but simply consider 
the condition that all forms of energy density must be non-negligible for an 
extended period, which we very conservatively estimate to be t = jq(H 0 _1 ), 
where i is of order unity, or a few Gyr. 

3.1 Discussion 

In the models we start with the phantom held energy density being dominated 
by the density of the ordinary matter. That is, we begin the evolution close 
to the equilibrium point F. Eventually the energy density of the phantom 
held dynamics ’’switches on” (the precise details depending on the shape of 
the phantom held potential V{<p) 0). The model then evolves towards the 
equilibrium point R. 

Numerical experiments were performed for different initial conditions. 
The results presented in Figure 1 are typical. In Figure 1, f \ and Q p h are 
only non-negligible for a very brief period, fib changes from unity to zero, 
as it moves away from the equilibrium point F towards the sink R, at an 
exponential rate (similarly for Q p h). For all numerical values (and for all 
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values of k) we find that t is negligible, and the physical time constraint 
cannot be satisfied. We note from Figure 1 that during the epoch in which 
i\ and Vtph are non-negligible, the curvature is also non-negligible, which 
is not consistent with current observations. By choosing particular initial 
conditions (fine-tuning), we can arrange for the curvature to be small. But 
no amount of fine-tuning can allow for a significant t. Hence these models 
are not physical. We emphasise that this physical time constraint is a very 
severe constraint and can be used to rule out many cosmological models. 

The behaviour to the past (which does depend on the sign of k) was briefly 
discussed above. The numerics confirm that in all cases O 5 crosses unity at 
a finite time in the past, whence Q p h becomes negative and eventually the 
variables diverge. Presumably additional terms must be added to the models 
to avoid this type of unphysical behaviour. 
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